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ABSTRACT 

We have constructed two types of analytical models for an isothermal hla- 
mentary cloud supported mainly by magnetic tension. The hrst one describes an 
isolated cloud while the second considers hlamentary clouds spaced periodically. 

Both the models assume that the filamentary clouds are highly flattened. The 
former is proved to be the asymptotic limit of the latter in which each hlamen¬ 
tary cloud is much thinner than the distance to the neighboring hlaments. We 
models reproduce main features of the 2D equilibrium model of 
for hlamentary cloud threaded by perpendicular magnetic held. 

It is also shown that the critical mass to hux ratio is M/<f) = {2tt^/G)~^, where 
M, $ and G denote the cloud mass, the total magnetic hux of the cloud, and 
the gravitational constant, respectively. This upper bound coincides with that 
for an axisymmetric cloud supported by poloidal magnetic helds. We applied the 
variational principle for studying the Jeans instability of the hrst model. Our 
model cloud is unstable against fragmentation as well as the hlamentary clouds 


show that these 
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threaded by longitudinal magnetic field. The fastest growing mode has a wave¬ 
length several times longer than the cloud diameter. The second model describes 
quasi-static evolution of filamentary molecular cloud by ambipolar diffusion. 

Subject headings: ISM: clouds; ISM: magnetic fields; magnetohydrodynamics 
(MHD); stars: formation 


1. Introduction 


Interstellar magnetic field plays significant roles on the dynamics of molecular clouds. 
Magnetic field supports molecular cloud against gravitational collapse if it is strong enough. 
It is well known that the mass to flux ratio, i.e., the amount of gas contained in a mag- 
netic flux tube, is the key parameter for the magnetic support. iMestel fc Spitzerl fll956l) 
pointed out from the viral analysis that magnetic force reduces gravity by a constant factor 
independent of the size of the cloud. They evaluated the critica l magnetic flux to support 
a cloud of mass, M, to be $c = tt^/GM ( Mestel fc ParisI 11984 1. The virial theorem pro¬ 
vides the correct scaling law but not quantitatively corr ect value for t he crit ical mass. One 
needs equilibrium models to evaluate the critical value. iMouschoviasI fll976alJbl) for the first 
time obtained a magnetohydrostatic equilibrium for a disk like cloud supported in part by 
magnetic field. He evaluated magnetic pressure and tension as well as self-gravity and gas 
pressure consistently. He confirmed that the mass to flux ratio is one of the key parameters. 
He also noticed that the equilibrium depends also on ratio of the magnetic pressure to the 
gas pressure and the flux function, i.e ., the amount of gas contained in each magnetic flux. 
Tomisaka. Ikeuchi fc Nakamura! (1l988l ) obtained a fitting formula for the critical equilibrium 
from their numerical results. It reads = (0.17)“^M\/G ~ 2tt\/GM. Finite electrical 
conductivity changes the flux function through ambipolar diffusion and induces quasi-static 
contraction. When the magnetic field is weaker than the critical one, the gravitational col¬ 
lapse continues to form stars. Thus magnetohydrostatic equilibria help our understanding 
of the gravitational collapse to form stars. 

Magnetohydrostatic equilibria depend on the configuration. We need to consider fila¬ 
mentary clouds in magnetohydrostatic equilibria since observed molecular clouds are often 
filamentary and associated with either longitudinal or perpendicular magnetic field. The 
magnetic field is parallel to elongation of the molecular c loud in the Ophiuchns star forming 
region while it is perpendicu lar in the Taurus ( see, e.g..lMoneti et al.]l9^ : Goodman et al 
1990 : Palmeirim et al. 2013 ) and in the Musca ( Perevra fc MagahaesI 20041 ). 


Stodolkiewczl fjl963|) and lOstrikerl (119641) obtained equilibrium model having the density 






































distribution, 


p{r) 

- '’“P + Si/O ■ 

(1) 
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for an isothermal filamentary cloud in the cylindrical coordinate, (r,ip,z), where Cg and G 
denote the isothermal sound speed and gravitational constant, respectively. Interestingly, 
the line density defined by 


A 


2nrp{r)dr = 



( 3 ) 


does not depend on the filament width, H. Thus, this model is neutrally stable against 
radial contraction. 


Stodolkiewca fjl963h showed that filamentary clouds having a larger line density can 
be sustained by purely longitudinal or partially helical magnetic field. These models are 
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1993 

. The fragmentation and its further evolution were studied extensivelv 

by 

Tomisaka 

1995) 

. Nakamura et al. 

(1995 

) and others. 


On the contrary, filamentary clouds permeated by magnetic held per pendicular to the 
axis have been studied little. Its structure was obtained only recently by iTomisakal (120141) 
referred to Paper I in the following. This is mainly because the structure is obtained only 
by extensive numerical computation solving the Poisson and Grad-Shafranov equations si¬ 
multaneously. The latter is highly nonlinear since it describes force balance in the directions 
parallel and perpendicular to the magnetic held. It requires high spatial resolution especially 
when the hlamentary cloud is condensed. Paper I showed that the maximum line density 
increases in proportion to the magnetic hux. In this paper we show analytical models which 
reproduce the main features of the 2D models obtained in Paper I. 


When the line density is higher than 2c^/G and thus supported mainly by magnetic 
held, the hlamentary cloud is battened like Italian pasta fettucine. The magnetic force is 
dominated by the magnetic tension. Then we can use the thin disk approximation to analyze 
the structure and stability of hlamentary cloud permeated by perpendicular magnetic held. 


This paper is organized as follows. In §2 we derive an analytical model for an isolated 
hlamentary cloud. The mass to hux ratio is uniform in the model. In §3 we use the Fourier 
series for hlamentary clouds arranged periodically. It is shown that the solution approaches 
to that obtained in 52 when the hlament thickness is much smaller than the distance to the 
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neighboring filaments. It is also shown that the ratio of the thickness to the distance can 
be regarded as the model parameter specifying the quasi-static evolution. In §4 we confirm 
that our ID models reproduce the main features of the 2D model obtained in Paper I. In §5 
the model obtained in §2 is proved to be unstable against fragmentation. In §6 we discuss 
the stability of the models obtained in §3. 


2. Equilibrium Model 


In this paper, we consider an isothermal flattened filamentary cloud having infinite 
length and infinitesimal thickness. This simplification is justified for observed filamentary 
clouds since their lengths are much larger than the widths and the magnetic fields are often 
perpendicular to cloud elongation on the sky. It is well known that self-gravitating clouds 
are flattened in the direction perpendicular to the global magnetic field. In the following 
we use the Cartesian coordinates in which the cloud is confined in the plane of z = 0 and 
elongated in the ^/-direction. Both the gas and magnetic fields are uniform in the ^/-direction. 

From the above assumptions we can express the gas density distribution as 


p = T.{x)5{z), 


(4) 


where D(a;) and 6{z) denote the surface density and the Dirac’s delta function, respectively. 
The Poisson equation reduces to 


dx'^ dz"^) 


0 


4:7iGTi{x)S{z). 


(5) 


We obtain the boundary condition. 


d(j) 


dz 


z=+e 



27iGT,{x). 


( 6 ) 


by integrating E quation fl5l) in the vertica l direction over the infinitesimal interval around 


z = 0 (see, e.g., iBinnev Sz Tremainell2008l. for the derivation). The symbol, s, denotes an 


infinitesimally small positive quantity. The gravity is given by 


9x 

9z 


d(t) 

dx ’ 
d(j) 

dz' 


(7) 

( 8 ) 


We assume for simplicity that the magnetic field has only the x- and ; 2 -components 
in the equilibrium. The condition for magneto-hydrostatic equilibrium requires that the 
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magnetic field should be force free and hence current free outside the disk, 

dB^ A-wJy^x) 

since both the gravity and pressure thereof are negligibly small. The symbols, c and Jy{x), 
denote the speed of light and the electric surface current running in the ?/-direction inside 
the cloud per unit length in the x-direction, respectively. The magnetic held is kinked on 
the disk, 

. N . N 27iJJx) , , 

B^{x,z = +e) = -B^{x,z = -e) = -. (10) 

c 


Note that the gravity and magnetic held can be perfectly aligned outside the cloud. 
Both of them are divergence and rotation free, i.e., V-^ = 0, Vx^ = 0, V-S = 0, and 
V X .B = 0, outside the cloud, z ^ 0. In the following we focus on the case that the magnetic 
held is proportional to the gravity. 


in 2 ; > 0 and 


Bx = 

a 

Vg' 

B, = 

a 

Vg' 

Bx = 

a 

B, = 

a 


9x1 

9zi 


(11) 

( 12 ) 

(13) 

(14) 


in z < 0, where a denotes a non-dimensional constant. Thi s assumption is equivalent to the 
constant mass to hux ratio (isopedic in the terminology of IShu fc Lilll997|) . 


B^(x,z = 0) 

“ = 

thanks to Equation ([6]). Hence, the cloud is subcritical (supercritical) when a > 1 (a < 1). 
When the mass to hux ratio is constant, the magnetic tension working on the cloud is 
proportional to the gravity. 




B^By 
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in the cloud. 


(16) 
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Equations ffT7|) simplifies the equation for magnetohydrostatic balance in the x-direction, 

d ( 


d r 2 \ d [Bl + Bl 
& (“‘d + Ts; 


dx 


Stt 


dz V An 


9xPi 


(18) 


where Cs denotes the isothermal sound speed. Substituting Equation (ITTll into Equation 
integrated over the inhnitesimal interval around = 0, we obtain 


(l - CK^) gjl + c; 


dT. 

dx 


= 0 . 


(19) 


The magnetic held reduces the effective gravity by a factor (1 — a^). See lShu fc Lil (119971) for 
the validity of this approximation. Equation (IT^ indicates that the cloud is in equilibrium 
only when a < 1 (supercritical). 


In the following we consider the Lorentz prohle, 

A a 


E = 


n x"^ At a? 


( 20 ) 


as a model for the battened hlamentary cloud. The line density and FWHM of this model 
cloud are A and 2a, respectively. Using the method of images we obtain the corresponding 
gravity, 


9 = 

The pressure force is given by 


f 9x) 


z 2GX / 


\ 9^ J 


\z\ + ( z + ay \ 

\z\ 


X 


( 21 ) 


.dS. 


2ac?Ax 


"dx 7r[x2 + (|^|+a)2]2' 

We obtain the condition for the magnetohydrostatic equilibrium, 

G\ (1 - a^) - cl = 0, 

by substituting Equations fl2UD through (12^ to Equation f[T2D . Equation 

A = 




2G y \2ny/GJ 


2G 


( 22 ) 

(23) 

is rewritten as 

(24) 


by the help of Equation flT5l) . When Equation (l2Tll holds, the model cloud settles in an 
equilibrium for any a. In other words, the width of the cloud is not specified by A or 4). This 
character is common for isothermal hlamentary cloud in equilibrium. 
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The magnetic field is expressed as 


B = 



$ 

Tr[x‘^ + (|2;| + a)2] 



(25) 


The magnetic field lines are straight and kinked on the plane of the flattened cloud. They 
look as if they emanate radially from the lines of {x,z) = (0, ±a). The strength decreases 
inversely proportional to the distance from the lines. 


We can estimate the density by assuming the hydrostatic balance in the vertical direc¬ 


tion, 


p{x,z) = p(x,0)exp{c/[0(a;,O)-0(x,z)]} , 
cf) = GX In -f- 2 :^ -f a^) . 


(26) 

(27) 


Here the gravitational potential, 0, is evaluated by the method of image and the magnetic 
force is neglected in the hydrostatic balance. The magnetic force, j x B, vanishes outside the 
disk since the electric current is confined in the disk. The magneti c force is perpe ndicular 
to the field inside the disk. Hence it compresses the gas disk (see IShu fc Lil Il997l. for the 
evaluation of the magnetic force acting on the vertical structure). The condition, 

/• + CXD 


E(U = 


p{x, z)dz, 


(28) 


gives us 


p{x,0) = 


T{GXc 


s^) 


V^^T^V^T (GAC72-I/2)’ 
where T denotes the gamma function. Accordingly we have 

Aa r(GAc72) (a;2 + 


p{x,z) = 


TT 0Fr ( GXc -^ - 1/2) (3,2 + ^2 + ^2) 


G\c: 


(29) 


(30) 


Figure [U denotes the density distribution and the magnetic field for A = 5, 10, and 20 
cl/G. Here the abscissa and ordinate are denoted in unit of a while the density is shown in 
unit of c\/{Ga). The cloud is more flattened for a larger A. 


3. Periodically Arranged Filamentary Clonds 

In this section we relax the uniform mass to flux ratio, H/B = const., assumed in the 
previous section. The mass to flux ratio should increase at the cloud center through the 
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Fig. 1.— The cross sections show the density by color and the magnetic field by lines. The 
npper left and right panels are for A = 5 and lOc^/G, respectively. The bottom one is for 
A = 20c^JG. 
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ambipolar diffusion if the cloud is supported in part by the magnetic held against gravity. 
On the contrary, the mass to hux ratio decreases outside the cloud. In order to take account 
of magnetic held outside clouds, we consider the situation in which hlamentary clouds are ar¬ 
ranged periodically on a plane. Again each hlamentary cloud is assumed to be geometrically 
thin. Then the surface density is given by the Fourier series, 

^ n=0 k / 

where ^ denotes the spacing between the hlaments. The 0th Fourier series coefficient, Oq, 
denotes line density per each hlament, A. We obtain the gravitational potential, 

27rGaO| , ^ Gan f 2nnx\ f 2nn \ 

(t){x,z) = —^|^|-^^cos(^-^Jexp(^- —\z\J (32) 


by solving the Poisson equation. The gravity is expressed as 


9x 


9z = 


2ttG / 2'!inx \ 

■—rl^an^in ^J exp 

f2nnx\ 


n=l 


27171 , 


2nG z 

■“TR 


E 

n=0 


a„ cos 


V ^ ; 


exp 


27771. 


Similarly the force free magnetic held is expressed as 


B = — 

" i\z\ 


bn sin i 

n=l ^ 


. f 27inx\ 

T~ J 


exp 


27777 , 


R 1 '^, f27777x\ 


exp 


27777 , 


(33) 

(34) 


(35) 

(36) 


The 0th Fourier series coefficient, bo, denotes the magnetic hux permeating each hlament, <F 
in the model for an isolated hlamentary cloud. 


The condition for the force balance. 


2^ _ 

dx 277 


fi'xE, 


(37) 


z=+e 


is rewritten as 


cRa, 


n 

m=0 


'mbn—m 


fE 

m=0 


^n—rrf 


(38) 
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This equation can be solved successively if once a = ( 01 , 02 ,...) and 60 are given. The 
Fourier series coefficient, is given by 




SttG G 


E 

m=0 


n—1 

^m^n—ra C^TIQ^ ^ ^ 


SttG 


n—m I 5 


(39) 


m=l 


if Oo, ..., o„, bo, ..., &n-i are known. In other words we can solve Equation fl371) if the mass 
distribution and the magnetic flux are given. 


In the following we consider the case. 


Oj — 


A (J = 0) 

2Aexp(—toj) (otherwise) ’ 


where w denotes the ratio of the hlament width to the interval. Then Equation 
to 


(40) 

reduces 


1 


E 

m=0 


bmbn—m 


2nXe-^^ {GX - cl) . 


(41) 


The right hand side of this equation indicates that the gravity overcomes the pressure force 
when A > cl/G. 


After some algebra we hnd a solution. 


$2 

47r^ 


r <F (j = o) 

\ 24) exp (—toj) (otherwise) ’ 
GX^ - clX. 


(42) 

(43) 


The mass to flux ratio is constant in this solution. The surface density and magnetic fields 
are expressed as 


E(a;, ^ = 0) 


B^{x,z = +£) 

B^{x,z = 0) 


A sinh w 


cosh w — cos 
271X 


$ sin 



cosh w — cos 
$ sinh w 


cosh w — cos 



(44) 


(45) 

(46) 
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respectively. See Appendix A for more details on the derivation. This solntion of w = ‘I'Kajl 
approaches to that given in §2 for a given a in the limit of £ —)■ cxd and hence tc —)■ 0. 

The solntion obtained in the previons paragraph is a critical one. When Bq is larger 
than the critical valne, 


bo > 2TTy/G 


ao 



(47) 


we obtain solntions in which the mass to flnx ratio is not nniform. Eqnation (H711 can be 
rewritten as 






(48) 


the right hand side of which denotes the maximum line density supported by magnetic flux, 
bo- The hlamentary clouds are conhned by the magnetic pressure of the neighboring clouds. 
Remember that the filamentary cloud is confined by uniform magnetic field in the 2D model 
of Paper I. Thus the model shown in this section is closer to the 2D model than that in the 
previous section. 

Figure [2] shows the magnetic field for A = 5c^/G and bo = 107rcf/\/G. The color denotes 
the density evaluated to be 


pix,z) 


7rGS(x)2 

2c2 



7rG'E{x)z 



(49) 


in the logarithmic scale. Equation denotes the equilibrium density distribution for an 
isothermal plane parallel disk having the surface density, S(x). The magnetic field is vertical 
to the disk plane near |a;| = £/2. 


Figure [3] shows the distribution of for models having various magnetic fluxes, = 
bo = 9, 10, 12, 14, 16, 18, and 20 ttcI/\/G. All the models have the same line density 
(A = bcl/G) and accordingly the same surface density distribution shown by the dashed 
curve. The value of w is fixed at w = 0.1. The magnetic field and the surface density 
profiles are very similar near the filament axis {x ~ 0). In other words almost the same 
amount of magnetic flux is confined in the filament and the mass to flux ratio is constant, 
Bz = 27r\/GE. When the magnetic flux is larger, the magnetic field is stronger outside the 
hlamentary cloud. When the magnetic hux is less than a critical value (bo < 97rc^i\^ in 
case of A = 5cf/G, see Eq. (1471) 1. we cannot construct an equilibrium model. When it is 
close to the critical value, the mass to hux ratio is nearly constant also outside the cloud as 
in the model shown in 52. 









0.4 



- 0.4 - 0.2 0.0 

xll 


0.2 0.4 


Fig. 2.— The cross section shows the logarithm of the density by color and the magnetic 
field by lines for A = 5c^/G and <h = lOncl/y/G. respectively. 



Fig. 3.— The solid curves denote the distribution of for A = 5cj./G. The surface density 
is denoted by the dashed curve. 
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Next we compare the models having the same A and $ but different w. Decrease 
in w mimics quasi-static contraction since the surface density at the cloud center increases. 
Figure m denotes the mass to flux ratio at the cloud center, 27r\/GSc/i?2,c; as a function of Sc, 
where Sc and B^ c denote the surface density and magnetic held at a; = 0, respectively. The 
ordinate is denoted in unit of c^/(Gt'). The red curve denotes the models having A = 6c^/G 
and <F = 407rc^/\/G while the black one does those A = 4c^/G and <F = 407rc^/\/G. As 
the surface density at the cloud center increases, the mass to hux ratio there increases 
monotonically and is saturated at 


27r\/GS 


S—^oo 



(50) 


Note that the mass to hux ratio tends to A/<F in the limit of tc —)■ cxo since both the surface 
density and magnetic held are uniform in the limit. Interestingly, however, the mass to hux 
ratio at the cloud center depends little on A for an intermediate w. 


When the line density is smaller than the critical value, c^g/G, the hlamentary cloud is 
conhned not by the self-gravity but by the magnetic held. This is because the gas pressure 
dominates over the gravity. We can construct such a model by assuming a small A and a 
large magnetic hux. They are similar to models with small Rq in paper I. 


4. Comparison with 2D model 

The thin disk model shown in §2 and 3 reproduces main features of the 2D model shown 
in Paper I. 

Equation fl2^ is essentially the same as Equation (45) of Paper I, 

= 2cl\ - GX^ + (51) 

OTT 

where R and ps denote the radius of the hlament and the gas pressure at cloud surface, 
respectively. The symbol, <hci, denotes a half of the magnetic hux permeating the hlament 
per unit length. Equation (15T|) was derived from the virial analysis and the left hand side 
denotes the pressure force acting on the cloud surface, which vanishes in our ID model. 
Diherences in the numerical factors come from the assumptions used. In this paper we 
neglected the hnite thickness of the cloud and hence pressure force acting in the ^-direction 
for simplicity. Thus the critical line density is evaluated to be Xcr = c^/G in case of no 
magnetic held (<F = 0) in this paper while it should be 2cl/G in paper I. Another diherence 
comes from the assumed mass to hux ratio; it is uniform in the model shown in §2 while it 
is not in Paper I. 
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The following equation, 


A = 


eg + (c! + ^4)^8) 
G 


1/2 


was derived from Equation (45) of Paper I as well as Equation 
). Equation (l5^ has an asymptotic form, 

^2 ^ _ 


A = ^ + ^ 


which resembles Equation (38) of Paper I, 


$ 1 

A^ax ^ 0-24^ + 1.66^. 


(52) 

is derived from Equation 

(53) 

(54) 


Equation (154]) is useful since it provides an upper limit on the line density supported against 
gravity by magnetic held and gas pressure. The upper limit increases in proportion to the 
magnetic hux contained in the cloud. Thus we examine the coefficient quantitatively. 

In paper I the mass to hux ratio is assumed to be the same as that of the hlament having 
uniform density, Pq, and threaded by uniform magnetic held, Bq. Then the line density and 
magnetic hux are denoted by 


A = TipoRl, 

*hcloud ~ 25oi?Q, 


(55) 

(56) 


respectively, where Rq denotes the radius of the uniform hlament. Note that the magnetic 
held surrounding but not permeating the hlament is not in the count of <hcioud- It should be 
also noted that <hcioud is twice as large as <hci used in Paper I. Hence, the diherential mass 
to hux ratio is expressed as 


dX 

= 


2A 


vrd), 


cloud 


1 


24)' 


4) 


cloud 


1/2 


Bz{x')dx'. 


(57) 

(58) 


Here the symbols, Bz{x) and dX/d^, denote the magnetic held on the plane of z = 0 and 
mass per unit magnetic hux, respectively. 

In order to evaluate the ehect of the mass to hux ratio distribution, we modihed the 
cross section of the initial hlament to be 

2IM 

/ u \ 

+ 


X 

Ro 


JL 

Ro 


2 


1 , 


(59) 
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where J\f > 0. The initial magnetic held is again uniform at Bo- Then the mass to hux ratio 
is expressed as 

dX _ r(A ^/2 + 3/2) A 
{J\f/2 + 1) ^ 

cloud 

where T denotes the gamma function. The mass to hux ratio is uniform in the cloud when 
A/" = 0. On the other hand, most of the gas is hlled in the magnetic held running at the 
center in case J\f = oo. Thus the parameter, A/", specihes the mass to hux ratio distribution 
as well as w in the model shown in §3. See Table [1] for the ratio of the mass to hux ratio at 
the center to the average. 

Figure [5] shows the mass to hux ratio at the cloud center, 2Tr'/GdX/d^ as a function 
of the density at the cloud center. Each curve denotes the locus of models having the same 
A and (3o, where (3o = PscI{Bq/S7i)~^ denotes the initial plasma beta in the ambient gas 
surrounding the hlamentary cloud. As in Paper I, the model clouds are assumed to be 
conhned by a very tenuous gas of which pressure is PsC^. The density at the cloud center is 
measured in unit of Bq/^Stic^). Thus the abscissa denotes the ratio of the gas pressure at 
the cloud center to the magnetic pressure in the region very far from the cloud, since the 
magnetic held is nearly uniform at the initial value, Bq, in our 2D model. The line density 
is specihed in unit of cl/G on the diagram. 

The index. A/", increases from A/" = 0.1 (left) to 10 (right) on each locus. This means 
that the hlamentary cloud is more centrally condensed and the central mass to hux ratio is 
higher when A/" is larger. The mass to hux ratio depends only a little on A in the models 
of /do = 0.01. This is because the gas is relatively cold and accordingly the gas pressure has 
only minor contribution to the cloud support. Note that the mass to hux ratio is slightly 
larger for a smaller A when pc is given. The maximum non-dimensional mass to hux ratio is 


2 $' 


<F 


cloud 


j\ j ^ 


(60) 


Table 1: The ratio of the central mass to hux ratio to the average. 


M 

(dA/d$),/(A/<Fcioud) 

0.0 

1.0000 

0.1 

1.0303 

0.2 

1.0598 

1.0 

1.2732 

2.0 

1.5000 

4.0 

1.8750 

10.0 

2.4610 
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Fig. 4.— The mass to flux ratio at the cloud center is shown as a function of the surface 
density at the cloud center for given A and <F, which are taken to be A = 4c^/G and 
$ = 87rc^/\/G for the black curve and A = 6c^/G and <F = 87rc^/\/G for the red curve, 
respectively. 



0.01 0.10 1.00 10.00 100.00 


0.01 0.10 1.00 10.00 100.00 


Fig. 5.— The mass to flux ratio at the cloud center is shown as a function of the density 
thereof in 2D models. The curves denote the loci of constant A of which values are labelled. 
The central density is denoted in unit of Bq/SttcI- The index, J\f, increases from left to right 
along the loci. The left and right panels denote the models of (3o = 0.01 and 0.1, respectively. 
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only slightly larger than unity when A > Sc^/G and /So = 0.01. 

When /So = 0.1, the gas is relatively warm and hence the pressure has a signihcant 
contribution to the cloud support. The non-dimensional mass to flux ratio is signiflcantly 
higher than unity especially when pc > Bq /although it is still less than three. This 
high mass to flux ratio is realized only when the cloud is relatively warm and the line density 
is relatively small. 

The maximum line density depends also on W, i.e., Amax = Amax(*^*ch Equation 

(46) of Paper is that for M = 1. When <l)ci and Cg are given, the maximum line density is 
lower for a larger A/". Increase in M mimics the quasi-static evolution of a molecular cloud 
by the ambipolar diffusion. As M approaches to zero, the maximum line density tends to 
the value given by Equation 


5. Stability Against Fragmentation 


We consider the stability of the isolated filamentary cloud against fragmentation, i.e., 
a sinusoidal perturbation in the 2 ;-direction. The gas is assumed to be still confined in the 
plane of ^ = 0 and we use the method of images to obtain the changes in the gravity and 
density consistently. The image density is assumed to have the form. 


gix,y,z) 

go{x,z) 

gi{x,y,z) 


go{x,z) + gi{x,y,z), 

2X6{x)6{z + a), 

6X{b) cos ky6{x)6{z + b)db, 


(61) 

(62) 

(63) 


Jo 

for calculating the change in the gravity in the upper half space of z > 0. The change in the 
gravity is evaluated as 


6g^{x,y,z =+e) = -G 


r’+oo 


(z - z')gi{x',y',z') 


= —Gb cosky 


[(x — x'y + {y — y'Y + 

5X{b) cos ks 


dx'dy'dz' 


'•+00 poo 


= —2Gkb cosky 


'o (a;2-I--f- 
6X{b)K^ {kVx^ + 62) 


dsdb 


db, 


(64) 

(65) 

( 66 ) 


\/a ;2 -f 

where Ki denote s the modified Bessel functio n of the 1st order. Here we used the mathe¬ 
matical formula (lAbramowicz fc Stegunlll965l ). 


KJz) = T Z/+- 


lA {2zy 


cost 


TT Jo (t2 + ;^2)^+1/2 


dt. 


(67) 
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Similarly we obtain 


5gx = —2Gk cos ky 


x5\{h)K^ {kVx^ + 62) 


6gy = 2Gk sin ky I 6X{b)Ko (^ky/x"^ + db. 


db, 


'0 

^OO 


The change in the surface density is expressed as 
SS(x,y,z = 0) = k_b^ r 


TT 


+ 62 


( 68 ) 

(69) 


(70) 


The magnetic held is assumed to be aligned with the gravity also in the perturbed 
state. This assumption is reasonable since the gas density is quite low and hence the inertia 
is negligibly small outside the cloud {z ^ 0). The magnetic held is still current free and 
does not extract any momentum from the cloud. Then the change in the magnetic held are 
expressed as 


= 


6By = 


6B, = 


a 


Vg 

a 

a 

77 


^Qxi 

^9yi 

Sgz, 


(71) 

(72) 

(73) 


in 2 : > 0. Both the gravity and magnetic held change with the time in proportion to 
exp{—iut), where u denotes the angular frequency of the perturbation. 


We use the displacement vector, 


^ = 


cos ky 
G sin ky 


exp(—icat), 


(74) 


for our stability analysis. The y-dependence of the displacement vector is chosen so that the 
the equation of mass conservation be the ordinary diherential equation with respect to x, 


^ ('CxSo) + k^yTlQ — 0, 


(75) 


where Sq denotes the surface density in the equilibrium. Also the equation of motion reduces 
to the ordinary diherential equation and algebraic equation, 


u;^j:o^x + cl^dJ:-{l-a^){j:o6gx + gx5J:) = 0 , 

—— CjME + (1 — Eo^Pj, — 0. 


(76) 

(77) 
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It is difficult to solve Equations fl75|l through (1771) simultaneously since they are as¬ 
sociated with Equation fl66|l . fl68|) and fl69|) . Instead of solving the equations we use the 
variational principle to evaluate First we rewrite Equations flTH]) and fl77|l into 


+ = 0, (78) 

—+ elk— -I- (1 — a^)5gy = 0, (79) 

by using the condition for magnetohydrostatic equilibrium, 

- (l - ^oQx = 0. (80) 

Note that the displacement is irrotational, 

+ k^x = 0, (81) 

ax 

since angular momentum extraction by magnetic held is not taken into account. We obtain 

— ca^Eo + ^y) + el-^ -(l — Eq + iy^9y) + ^ (csCx<^E) = 0, (82) 

by taking the sum of the products of c^^E/Eq and Equation (1751) . ^a;Eo and Equation fl78|) . 
^yEo and Equation d79|) . We obtain 


u 


2 



(83) 


(84) 


by integrating Equation fl8^ and submitting = 0 at x = ±cxo. The right hand side of 
Equation flHT)) gives a lower bound for the eigenvalue, when an arbitrary perturbation is 
substituted. If it is negative for a given perturbation, the hlament is unstable. 

We evaluate the right hand side of Equation flM)) as a functional of the change in the 
image density, 6X{b). The change in the gravity is evaluated by numerical integral of fl68|) . 
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and (EH]). The change in the snrface density is evalnated by the nnmerical integral of flTOjl . 
The displacement, is obtained by solving the ordinary differential eqnation, 


d d 

dx V ^0 / dx 




k So^a; — 0, 


(85) 


which is obtained by combining Equations fl75|l and flHTl) . The displacement, is obtained 
simultaneously when we solve Equation ([85]). The above mentioned procedure allows us to 
evaluate the right hand side of Equation fIMD as a functional of 5\. See Appendix B for more 
details. 


As shown in the previous section, FWHM of our model cloud (=a) is not specihed by 
the condition for magnetohydrostatic equilibrium. If the hlamentary cloud has proper line 
density and magnetic flux, it can be settled in equilibrium at any width. This means that our 
model cloud is neutrally stable for perturbation having A; = 0 . As in the case of longitudinal 
or helical magnetic held, we use the non-dimensional wavenumber, /ca, in our analysis. 

In this paper we consider two types of trial functions for the variational principle. Type 
I trial function is expressed as 


5\i{h) = 5{h-l3a), ( 86 ) 

where (3 is chosen to minimize the value of —Type II is expressed as 

N 

5X^{b) = Idjo), (87) 

where a set of 9j is chosen to minimize the value of for hxed (dj. 

The left panel of Figure El shows the growth rate obtained by applying type I trial 
function. The abscissa is the wavenumber in unit of a~^ while the ordinate is the growth 
rate in unit of [cs/aY- The right panel of Figure El shows that obtained by applying type 
II trial function in which the values of [3i are taken to be /3 = (1.0,1.4,1.8, 2.2, 2.6). Our 
model hlamentary clouds are unstable against fragmentation as well as the hlamentary clouds 
threaded by longitudinal magnetic held. The most unstable mode has the wave number, 
fcmax — 0.8a“^, and the growth rate, |a;| ~ O.bc^a”^, irrespectively of A. 


6. Discussions and Summary 

As shown in the previous sections, our ID model based on the thin disk approximation 
reproduces the main features of 2D model shown in paper I. The maximum mass sustained 
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by magnetic field is roughly proportional to the magnetic fiux. The critical mass to fiux 
ratio depends on the mass loading, i.e., the distribution of T,/B. When the mass loading 
is uniform (S/i? = const.), it is slightly larger than X/B^ > The same critical 

ratio is obtained f o r disk s in equilibrium and slabs stable against fragmentation. Remember 
that IStrittmatterl fll966l ) obtained a similar value from the viral analysis on the magneto¬ 
hydrostatic equilibrium. The critical ratio seems to depend little on the cloud geometry. 
Remember that filamentary clouds with a longitudinal field are unstable against fragmen¬ 
tation even if the magnetic field is very strong. The instability is due to the fact that the 
mass to fiux ratio is infinitely large since the clouds are highly extended in the direction of 
magnetic field. 

Given the critical mass to fiux ratio depends little on the cloud geometry, then it can be 
applied to the stability of periodically arranged filamentary clouds. They are almost isolated 
each other in the limit of w —?• 0. Thus they are unstable against fragmentation as well as 
the isolated filamentary clouds since the mass to fiux ratio is larger than the critical value at 
the cloud center (Sc > 2tt'/GBz^c)- When w —>■ cx), our ID r nodel reduces to a slab of which 
stability was already investigated bv iNakano fc Nakamura! (119781 ). The slab is stable as far 
as S < 2tt\/GB. In short the models shown in Figured] are stable/unstable for small/large 
value of Sc. We surmise that the models of Sc < 2'Ky/GB,, c are stable against fragmentation. 
Most of the 2D models shown in Figure [5] are also likely to be stable against fragmentation 
since Sc < 27r\/GR,,c. 

The above argument brings us an interesting result. Filamentary clouds supported 
by longitudinal magnetic fields are unstable against fragmentation while those supported by 
perpendicular fields can be stable. The latter can collapse quasi-statically through ambipolar 
diffusion, while the former cannot. 

Our models are still idealistic since turbulence and other dynamical effects are not taken 
into account. Nevertheless, they provide physical insights on the dynamics of filamentary 
clouds. Filamentary clouds with perpendicular magnetic field are likely to be able to sustain 
their forms as far as they are subcritical at the cloud center. They may fragment after the 
mass to fiux ratio e xceeds the critical value by a mbipolar diffusion or by mass accretion along 
the magnetic field f|Heitsch &: HartmannI 1201411 . 
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A. Fourier Series 


We used the mathematical formula, 

CO 


a"" cos nx = 


n=0 


a cos X — a 
1 — 2a cos a: + ’ 


(Al) 


for |a| < 1 in order to derive Equation from Equations (l3TD and fHOj) . The gravitational 
potential, 0, and the ^/-component of the vector potential. Ay, are expressed as 


2nGX. 


z — 


^ GX /27inx\ 

— exp 

n=l ^ ' 


27r, 

-n\ w — —\z\ 


2 vr, 

-n [ w — —\z\ 


^x <h 27rnx 

Ay = ^+> —sm^—exp 

i Tin t 

n=l 

respectively, for the solution specihed by Equations (HT|) through 
pressed as 


= -GAlog 


27r|2:|\ f 271X 

cosh w H- - — I — cos 


Ay = — tan 

71 


-1 


{TlX 


tanh I 2 + — 


respectively, since 


a'" 1 / 2 \ 

cosnx = —-log (1 — 2acosa; + a j , 


E ll 

n ' 2 

n=l 

E d 

— sm nx = tan 

n 


n=l 


asinx 
1 — a cos X j ’ 


for |a| < 1. 


(A2) 
(A3) 

They are also ex- 
(A4) 

(AS) 


(A6) 

(A7) 


B. Perturbation Equations 

For a given we obtain and by the following procedure. 

First we rewrite Equations fISS]) and flSTD as 

dlnSo, ,, 5E 


_ 

dx 

^ = -kL- 

dx 


(Bl) 

(B2) 
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We consider the case in which and are symmetric and antisymmetric with respect 
to X, respectively. This choice is rational since the nnpertnrbed state is symmetric and we 
are interested in the gravitational instability. An eigenmode pertnrbation should be either 
symmetric or antisymmetric and the fragmentation of a hlamentary cloud is symmetric with 
respect to x. Thus we obtain the boundary condition, ^ 3 , = 0 at x = 0. We assume that the 
displacement diminishes in the form, 

^x oc exp {—kx ), (B3) 

oc exp {—kx ), (B4) 

in the limit of x —)■ +cxo. This choice is based on the fact that the relative change in the 
surface density is proportional to exp(—/cx) in the limit. Equations flB3p and flB4p are valid 
also when (5S/So decreases more steeply with increase in x. 

We integrate Equations fIBip and flB2p from x = 0 to Xout with stepsize Ax = 5 x 10“^ by 
the 4th order Runge-Kutta method to obtain two sets of solutions. One has initial condition 
(Cx) d^x/dx) = (0, 0). The other is the solution of Equations flBlD and (IB2p without the source 
term, i.e., the homogenous one and has the initial condition (^ 3 ,, d^x/dx) = (0,1). We obtain 
the solution satisfying the boundary conditions at x = 0 and Xout as a linear combination 
of them. The outer boundary is set at Xout = min(15/c“^,40). The y-component of the 
displacement, is derived from Equation flHip . 

Now we can evaluate the right hand side of Equation fIMp for a given 6X{b). Consequently 
the growth rate is evaluated to be a function of b in type I trial function. We searched for 
b which maximizes the growth rate in the interval of 0.55a < b < 2.5a with the interval 
Ab = 0.05a. Both the denominator and delimiter of the right hand side of Equation fl84p 
are the quadratic expressions of 6*1,..., 6*5, when type II trial function is used. We used 
LAPACK, a mathematical library for linear algebra installed in IDE, the interactive data 
language, for minimizing 

We used the public Fortran subroutines downloaded from Prof. Jian-ming Jin’s web page 
(jin.ece.illinois.edu) to obtain the numerical values of the modihed Bessel functions. 
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Fig. 6.— The growth rate is shown as a function of the longitudinal wavenumber. It is 
derived by applying the type I and type II trial functions to the variational principle, in left 
and right panels, respectively. See text for further details. 







